o 



u 
< 



PACS numbers: 



OU-TAP-268 
YITP-06-06 

Angular Trispectrum of CMB Temperature Anisotropy from Primordial 
Non-Gaussianity with the Full Radiation Transfer Function 

Noriyuki Kog 

Department of Earth and Space Science, Graduate School of Science, 
Osaka University, Toyonaka 560-0043, Japan and 
Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

Eiichiro KomatsrQ 

Department of Astronomy, The University of Texas at Austin, 1 University Station, C1400, Austin, TX 78712 

We calculate the cosmic microwave background (CMB) angular trispectrum, spherical harmonic 
transform of the four-point correlation function, from primordial non-Gaussianity in primordial 
curvature perturbations characterized by a constant non-linear coupling parameter, /nl- We fully 
take into account the effect of the radiation transfer function, and thus provide the most accurate 
estimate of the signal-to-noise ratio of the angular trispectrum of CMB temperature anisotropy. 
We find that the predicted signal-to-noise ratio of the trispectrum summed up to a given / is 
approximately a power-law, (S/N)(< I) ~ 2.2 x 1CP 9 /niJ 2 > u P to the maximum multipole that we 
have reached in our numerical calculation, I = 1200, assuming that the error is dominated by cosmic 
variance. Our results indicate that the signal-to-noise ratio of the temperature trispectrum exceeds 
that of the bispectrum at the critical multipole, l c ~ 1500 (50/|/nl|). Therefore, the trispectrum of 
the Planck data is more sensitive to primordial non-Gaussianity than the bispectrum for |/nl| > 50. 
' We also report the predicted constraints on the amplitude of trispectrum, which may be useful for 

other non-Gaussian models such as curvaton models. 
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\ Inflation has been the standard paradigm for the origin of cosmological fluctuations. While simple inflationary 
£vq • models based on a slowly-rolling scalar field are unable to generate the detectable level of primordial non-Gaussianity, 
a large class of models predict much stronger signals (see Q for a review). Non-Gaussianity of the primordial 
fluctuations thus plays an important role in testing and constraining inflationary models. As different statistical 
methods are sensitive to different aspects of non-Gaussianity, one should explore a variety of methods in order to 
maximize our sensitivity to primordial non-Gaussianity. Of which, the higher-order correlation functions such as the 
Q H l three- and four-point correlation functions, or their harmonic counterparts, the bispectrum @,0,0,IEIE@iHil3 
and trispectrum [Tol IllL Il2t Il3j | , have been actively investigated in the literature as a powerful probe of primordial 
■ non-Gaussianity. 

Compared with progress in theoretical calculations of the angular bispectrum of CMB temperature 0, H| and 
polarization [(| anisotropy, that for the angular trispectrum has been a little bit behind. The first calculation done by 
[l3| did not include the full effect of radiation transfer function caused by acoustic physics at the surface of last scatter. 
• i-h , In this paper, we calculate the angular trispectrum of CMB temperature anisotropy, fully taking into account the 
1 radiation transfer function. Our estimate of the signal-to-noise ratio of the trispectrum should thus improve accuracy 
of the previous estimate based on an approximate method. We then compare the signal-to-noise ratios of primordial 
■ - - > non-Gaussianity from the trispectrum and bispectrum. 

Let us briefly review the formalism of the CMB angular trispectrum, following [T2I . We decompose temperature 
anisotropy on the sky, ST/T, into spherical harmonic coefficients, a; m , as 

ST 

= ^2 a lmYl m (n). (1) 
Ira 

Statistical isotropy of the universe requires an n-point correlation function be rotationally invariant; thus, for n = 4 
one obtains 

(a hmi a l2m2 a hma a Umi ) = ^(-1) M ( £ ^ J ( ^ ^ ^ J l££(£), (2) 
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where T/^fX) is the angular averaged trispectrum, L is the length of a diagonal that forms triangles with li and 

12 and with l 3 and £4, and the matrix is the Wigner 3-j symbol, which guarantees that two sides and the diagonal 
form a triangle, \h — h\ < L < l\ + I2 and \h — h\ < L < I3 + 14. Parity invariance also requires h + I2 + L = even, 

13 + I4, + L = even, mi + — M = 0, and TO3 + + M = 0. These conditions determine the number of possible 
configurations. 

The trispectrum generically consists of the connected, T Cl and unconnected, Tg, part: 



Tl£(L)=Tj l ±(L)+T G \±(L), (3) 
The former contains non-Gaussian signatures, while the latter contains only the angular power spectrum, C'i = (|a; m | 2 ), 

T G \±{L) = (-l^+V^i + l)(2l 8 + i)c h c h s hl2 s hh s L0 



+ {2L + l)C h C l2 [(-iy i+h+L S lll3 5 hli +S hl J hh 
Using permutation symmetry, one may write the connected part of the trispectrum as 
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where 



P!±(L)=t\±(L) + (-lf 



•2L+h+l2+h+hj2l 



L+h+h+ht 



mm- 



(6) 



Here, the matrix is the Wigner 6-j symbol, and t ; ^(£) is called the reduced trispectrum, which contains all the 
physical information about non-Gaussian sources. 

We parameterize primordial Bardeen's curvature perturbation, during the matter-dominated era in the usual 
form as H [13 



$(*) = $ L (x) + /nl H(x) - (*l(x))] + h*l(x), 



(7) 



where $l is the linear Gaussian part and /nl and fi are the non- linear coupling parameters. (Note that $ is a 
perturbation in the (i, i)-component of the metric.) The current observational constraint on /nl is —58 < /nl < 134, 
which is from the analysis of the angular bispectrum of the WMAP data 0, 0] . No constraint on /nl from the 
angular trispectrum is currently available. It is therefore important to obtain an accurate estimate of the signal-to- 
noise ratio of the angular trispectrum of primordial non-Gaussianity expected from the WMAP as well as Planck 
experiments 0|. While we take /nl to be a constant independent of coordinates or wavenumbers, in general /nl 
depends on scales, and different inflationary models predict different dependence of /nl on wavenumbers [j. The 
post-inflationary evolution of $ due to second-order metric perturbations, which must exist in any models of the 
standard cosmology, also generates wavenumber-dependent /nl [!]• Nevertheless, a constant /nl model is still useful 
for estimating sensitivity of CMB experiments to the amplitude of non-Gaussianity. Measurements of non-Gaussian 
fluctuations are usually quite challenging as non-Gaussianity is very small, which makes detection of wavenumber- 
dependent features even more challenging. Since currently there is no prediction for sensitivity of CMB experiments 
to /nl from the angular trispectrum of CMB temperature anisotropy with the full radiation transfer function taken 
into account, we shall adopt a constant /nl to explore the signal-to-noise ratio of the trispectrum as a function of the 
maximum multipoles measured by observations. In the future one may extend our approach to include scale-dependent 
/nl, following the method given in Q, for instance. 

The primordial fluctuations yield temperature anisotropy as 



4tt(- 



d 3 k 
(2^)3 



^(k)g Tl (k)Y^(k), 



(8) 



where gri(k) is the radiation transfer function of adiabatic fluctuations, which can be calculated numerically by the 
CMBFAST code 0|- Using this relation, we obtain the formula of the reduced angular trispectrum as 



tL hlS L ) = / 7 i dr i r 2 dr 2 F L (r 1 ,r 2 )ai 1 (r 1 )f3i 2 (r 1 )a h (r 2 )(3u(r2)hi lL i 2 h hLU 
r 2 dr A 2 (r)/3; 4 (r) [/^ (r)$ 3 (r) + Ph(r)m 3 (r)] h hL i 2 h hLh , 



(9) 
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where 



F L {n,r 2 ) 


2 
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f k 2 dk P 9 (k)gri{k)ji(kr), 


IH(r) 


= *J 


[ k 2 dk f29Ti{k)ji(kr), 



(10) 

(11) 

(12) 
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and 



_ - / (2f 1 + l)(2i 2 + l)(2£ + l) ( h h L\ n ,, 

h hu 3 = y ^ J ■ (14) 

We have defined the power spectrum of the primordial curvature perturbation as 

($ L (fc)$£(fc')> = (2vr)V 3 )(fc - k')P*(h), (15) 

where P$ cx k n ~ A and n = 1 for a scale-invariant spectrum. 

The signal-to-noise ratio of the temperature trispectrum is given by 

Note that we have assumed that the error is solely due to the cosmic variance. When detector noise is included, Ci in 
the denominator should be replaced by C; + Ni/Wi, where 2Vj is the noisebias and Wi is the window function. We set 
the normalization of P$(fc) by the height of the first acoustic peak of the WMAP data; 1(1 + l)Ci/2ir = (74.7/zK) 2 at 
I = 220, and adopt the cosmological parameters of the concordance ACDM model with a scale-invariant primordial 
spectrum for h = 0.72, Q b = 0.047, Q A = 0.71, Q m = 0.29, and r = 0.17 While we shall evaluate (S/N) 2 
using the trispectrum that uses the radiation transfer function [equation @] below, let us first estimate an order- 
of-magnitude of (S/N) 2 and its dependence on the maximum multipole, l max , using the Sachs-Wolfe approximation 
valid at low multipoles (I <C 100), gn(k) ~ — j;(fcr*)/3, where r* is the comoving distance to the surface of last 
scatter. This form of gn(k) then gives 

t\±(L) « 9C™C™ [4/ 2 L Cf w + h (CfT + Cg")] h hL i 2 h h LU, (17) 

where 

C ? W = ^ J k 2 dkP <s> (k)j 2 (kr r ) = (18) 

is the angular power spectrum in the Sachs- Wolfe approximation for a scale-invariant primordial spectrum of and 
A ks 6 x 10~ 10 . For simplicity we take into account only the first terms in equations JS} and Moreover, let 
us consider only L = 1 modes because the collapsed configurations, which correspond to the low-L modes, are the 
dominant modes for /nl model [equation lj7j)]. This property also holds for the bispectrum For f 2 = we have 

t\±(L = 1) « WLClZcf^Cf^h^h^. (19) 

Since \h — 1| < 1% < l\ + 1 from the triangle condition and l\ > ?2j ^Zi,i.z 2 i s non-zero only for Zi — 1 = l 2 . Hence 
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FIG. 1: Signal-to-noise ratio squared, (S/N) 2 , of the angular trispectrum summed over all L, L < 10, and only L = 1 modes, 
respectively, as a function of the maximum multipole, i max . Here, L denotes the multipole of the diagonal of a trispectrum 
configuration. We assume /nl = 1 and /2 = 1 and use the full radiation transfer function. This figure shows that the summation 
over L needs to be done only up to L ~ 10. 



which gives 



N J ~ ^ 3C sw%w c sw c sw 



27V 1 ^ 



4! 3Cf w Cf w 

~ ^/nlCsxi (21) 

is the maximum multipole. This result indicates that (S/N) 2 of the trispectrum strongly depends on 
both / NL and Z max , (S/N) 2 ri cx (/NL^max) 4 , which is stronger than that of the bispectrum, (S/N)l { cx (/NL^max) 2 0- 
We have confirmed that ZJ^ax dependence still holds when all the terms (still in the Sachs- Wolfe approximation) are 
included. The constant of proportionality is however about 20 times larger than equation 1|21[1 . This is because there 
are other 11 terms due to the permutation symmetry as seen in equations JHJ and iffijl. and summing over L > 1 modes 
further increases (S/N) 2 by a factor of two. We shall show below that /^ax dependence also holds very well for the 
full calculations. 

As the trispectrum calculation involves the summation over five multipoles, computation takes too long to let us 
go beyond i ma x ~ 100; however, we find that almost all the contribution actually comes from L < 10 as shown in 
Fig. n where we use the full radiation transfer function. It is therefore sufficient to perform the summation over the 
diagonal, L, only up to 10, which makes the computational time scale as Z^ax instead of Z^axi gi vm g a huge saving in 
computational time. Figure |2] shows the predicted signal-to-noise ratio squared, (S/N) 2 , for the trispectrum (summed 
over L < 10) as well as the bispectrum as a function of the maximum multipole, ^ max , with the full radiation transfer 
function. We have assumed a fiducial value of /nl = 50 and ji — 1 in which case the effect of /2 is negligible. Our 
result suggests that the expected signal-to-noise ratio of the primordial trispectrum is approximately given by 



N 



5 x 10- 18 /nlC (22) 



or (S/N) tl i ~ 2.2 x 10 9 (/NiJmax) 2 - The minimum detectable /nl by the trispectrum at the 1-a level is therefore 
/nl ^ 85, 42, 21, 14, and 11 for 2 max = 250, 500, 1000, 1500, and 2000. This may be compared with the minimum 
detectable / NL by the bispectrum, / NL ^ 23, 12, 6, 4, and 3 for l max = 250, 500, 1000, 1500, and 2000. (Note again 
that we assume that the CMB bispectrum and trispectrum are cosmic-variance limited up to l max -) Here, we have 
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FIG. 2: Predicted signal-to-noise ratio squared, (S/N) 2 , of the angular trispectrum and bispectrum with the full radiation 
transfer function for /nl = 50 and fa = 1. Power-law fits are also shown. Note that the power-law fits break down at 
'max > 3000, where the gravitational lensing effects become important: (S/N) 2 ceases to grow at Z max ~ 3000 0-13 



used a power-law fit to (5/JV)bi — 1.7 x I0~ 4 l/NL^max- The power-law fits must break down at l max > 3000, where 
(S/N) 2 does not grow due to the gravitational lensing effects increasing the cosmic variance 

These results might seem to imply that the temperature trispectrum is always less sensitive to /nl than the 
bispectrum; however, as (S/N) tr [ depends on /nl more strongly than (S/N)^, the ratio also depends on /nl: 

(S/N) tli |/nl| I max /no\ 
{S/N) hi ~ 52 1500' 1 ' 

The trispectrum actually becomes more sensitive to primordial non-Gaussianity than the bispectrum for Planck-like 
experiments probing Z max > 1500 and |/nl| 50. Incidentally, the trispectrum cannot measure the sign of /nl? a $ 
it depends on / NL . While we have considered the trispectrum of temperature anisotropy only, the signal-to-noise 
ratio should increase when polarization is included in the analysis. The number of possible combinations of the 
temperature and _E-polarization field for the trispectrum is 5, while that for the bispectrum is 4; thus, one expects 
similar improvements in the signal-to-noise ratio for the bispectrum and trispectrum. The calculation including the 
temperature and polarization bispectra has shown that the signal- to- noise ratio increases roughly by a factor of 2 |(| . 
We expect a similar level of improvement for the trispectrum as well. 

Equation 11' ill) holds only for the particular non-Gaussian model given by equation (JJJ; however, the amplitude of 
the trispectrum may be related to that of the bispectrum in a very different way for other models. For example, the 
curvature perturbation may be given by 

$(*) = $l(x) + /, [r) 2 (x) - ( V 2 (x))] , (24) 

instead of equation J7J, where r\ is another fluctuating field that is totally uncorrelated with $l- This form of non- 
Gaussianity may arise from a particular configuration of curvaton models [T^. I20I ] . (The form of the CMB angular 
bispectrum arising from this model has been derived in Appendix C of [Toj^ In this model the ratio of the amplitude 
of trispectrum to bispectrum is very different from that for equation Q . Boubekeur and Lyth [21| proposed to use 
the following parameterization for the trispectrum amplitude of primordial curvature perturbation in the comoving 
gauge, C, which is related to Bardeen's curvature perturbation during the matter era as ( = (5/3)$: 

(C(fci)C(*a)C(*3)C(*4)> = (27r) 3 <5 (3) (fc 1 + k 2 + k 3 + fc 4 )T NL [P c {k 1 )P c (k 2 )P c (\k 1 + fc 4 |) + (11 distinct permutations)] . 

1-1 (25) 

Equation J7J) with f 2 = gives tnl = (6/nl/5) 2 , while equation (|2*4"|l gives different predictions |2l|. Our results 
therefore give the minimum detectable tnl at the 1-a level as tnl — 10400, 2500, 640, 280, and 170 for i max = 250, 
500, 1000, 1500, and 2000. These results should be useful for estimating sensitivity of CMB experiments to a variety 
of non-Gaussian models giving rise to non-zero trispectrum. 
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Our predictions show that the angular trispectrum measured from WMAP is unable to detect primordial non- 
Gaussianity parameterized by /nl as in equation J7J), given the current constraint on /nl from the bispectrum. 
However, the trispectrum is at least as powerful as the bispectrum for detecting primordial non-Gaussianity in the 
Planck data if |/nl| ~ 50, which is still allowed by the WMAP data. For other models, no detection of significant 
trispectrum on the WMAP (Z max ~ 250) and Planck data (Z max ~ 1500) would imply tnl £ 2 x 10 4 and 560 at the 
2-er level, respectively. 
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